AEROSPACE  REPORT  NO. 
TOR-0091  (6078)-2 


A241  853 


Relaxation  Processes  of  Vibrationally  Excited  H2O 
in  the  Mesosphere  and  Thermosphere 


Prepared  by 
V  I.  LANG 

Space  and  Environment  Technology  Center 
Technology  Operations 


1  September  1991 


•:  ?  r '  r 


nr; 


i  -J  v.1  i  j  ^ 


'Pi  A  ' 


Prepared  for 


PHILLIPS  LABORATORY 
GEOPHYSICS  DIRECTORATE 
Hanscom  AFB,  MA  01731 


SPACE  SYSTEMS  DIVISION 
AIR  FORCE  SYSTEMS  COMMAND 
Los  Angeles  Air  Force  Base 
P.  0.  Box  92960 
Los  Angeles,  CA  90009-2960 


Engineering  and  Technology  Group 


91-13300 


THE  AEROSPACE  CORPORATION 

El  Segundo,  California 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED 


TECHNOLOGY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  “architect-engineer”  for  national  security  programs, 
specializing  in  advanced  military  space  systems.  The  Corporation’s  Technology  Operations  supports  the 
effective  and  timely  development  and  operation  of  national  security  systems  through  scientific  research 
and  the  application  of  advanced  technology.  Vital  to  the  success  of  the  Corporation  is  the  technical  staff’s 
wide-ranging  expertise  and  its  ability  to  stay  abreast  of  new  technological  developments  and  program 
support  issues  associated  with  rapidly  evolving  space  systems.  Contributing  capabilities  are  provided  by 
these  individual  Technology  Centers: 

Electronics  Technology  Center:  Microelectronics,  solid-state  device  physics,  VLSI 
reliability,  compound  semiconductors,  radiation  hardening,  data  storage  technologies, 
infrared  detector  devices  and  testing;  electro-optics,  quantum  electronics,  solid-state 
lasers,  optical  propagation  and  communications;  cw  and  pulsed  chemical  laser 
development,  optical  resonators,  beam  control,  atmospheric  propagation,  and  laser 
effects  and  countermeasures;  atomic  frequency  standards,  applied  laser  spectroscopy, 
laser  chemistry,  laser  optoelectronics,  phase  conjugation  and  coherent  imaging,  solar 
cell  physics,  battery  electrochemistry,  battery  testing  and  evaluation. 

Mechanics  and  Materials  Technology  Center:  Evaluation  and  characterization  of  new 
materials:  metals,  alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of 
carbon;  development  and  analysis  of  thin  films  and  deposition  techniques; 
nondestructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  development  and  evaluation  of  hardened  components; 
analysis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures;  launch 
vehicle  and  reentry  fluid  mechanics,  heat  transfer  and  flight  dynamics;  chemical  and 
electric  propulsion;  spacecraft  structural  mechanics,  spacecraft  survivability  and 
vulnerability  assessment;  contamination,  thermal  and  structural  control;  high 
temperature  thcrmomcchanics,  gas  kinetics  and  radiation;  lubrication,  and  surface 
phenomena. 

Space  and  Environment  Technology  Center:  Magnetospheric,  auroral  and  cosmic  ray 
physics,  wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric  and 
ionospheric  physics,  density  and  composition  of  the  upper  atmosphere,  remote  sensing 
using  atmospheric  radiation;  solar  physics,  infrared  astronomy,  infrared  signature 
analysis;  effects  of  solar  activity,  magnetic  storms  and  nuclear  explosions  on  the  earth’s 
atmosphere,  ionosphere  and  magnetosphere;  effects  of  electromagnetic  and  particulate 
radiations  on  space  systems;  space  instrumentation;  propellant  chemistry,  chemical 
dynamics,  environmental  chemistry,  trace  detection;  atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and  radiative 
signatures  of  missile  plumes,  and  sensor  out-of-field-of-view  rejection. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


la.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 

2a7  SLCURITY  CLASSIFICATION  AUTHORITY " 


REPORT  DOCUMENTATION  PAGE 

1 1b  restrictive  MARKINGS 


2b.  DECLASSIFICATION/DOWNGRADING  SCHEDULE 
~4  "PERFORMING  ORGANIZATION  REPORT  N'JMBER(S) 

TOR-0091(6078)-2 


3.  DISTRIBUTION/AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  unlimited. 

5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


6a.  NAME  OF  PERFORMING  ORGANIZATION  1 6b.  OFFICE  SYMBOL  7a.  NAME  OF  MONITORING  ORGANIZATION 
,  .  (If  applicable) 

The  Aerospace  Corporation  Space  Systems  Division 

Technology  Operations _ _ . 

6o.  ADDRESS  (City,  State,  and  ZIP  Code)  7b.  ADDRESS  (City,  Stale,  and  ZIP  Code) 

El  Segundo,  CA  90245-4691  Los  Angeles  Air  Force  Base 

Los  Angeles,  CA  90009-2960 

'  oT^AmToF  FUNDING, SPONSORING  8b.  OFFICE  SYMBOL.  9,  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

OBG.TOTION  F04701-88-C-0089 

Phillips  Laboratory 

8c.  ADDRESS  (City  State,  and  ZIP  Code)  10.  SOURCE  OF  FUNDING  NUMBERS 

Geophysics  Directorate  program  [project  (task 

HanscomAFB,MA  01731  ELEUENT"°-  N°'  m 

11.  TITLE  (Include  Security  Classification) 

Relaxation  Processes  of  Vibrationally  Excited  H2O  in  the  Mesosphere  and  Thermosphere 


ro.  SOURCE  OF  FUNDING  NUMBERS 
PROGRAM  I  PROJECT  I  TASK 

ELEMENT  NO.  NO.  NO. 


WORK  UNIT 
ACCESSION  NO. 


12.  PERSONAL  AUTHOR(S) 

Lang,  Valerie  I. 


13a.  TYPE  OF  REPORT 


16.  SUPPLEMENTARY  NOTATION 


13b.  TIME  COVERED 
FROM _ TO 


14.  DATE  OF  REPORT  (Year,  Month,  Day)  15.  PAGE  COUNT 

1991  September  1  51 


COSATI  CODES 


GROUP  SUB-GROUP 


18.  SUBJECT  TERMS  f Continue  on  reverse  if  necessary  and  identify  by  block  number) 

H2O  (Water)  Upper  Atmosphere 

Non-Local  Thermodynamic  Equilibrium  Vibrational  Relaxation 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 


A  review  of  rate  constants  for  the  relaxation  of  vibrationally  excited  H2O  is  presented.  Tills  evaluation  is  the 
second  in  a  series  of  reports  intended  to  improve  the  state-of-know ledge  regarding  the  chemical  kinetics  of  the 
mesosphere  and  them,  sphere.  The  vibrational  levels  discussed  are  limited  to  those  presumably  populated  at 
altitudes  5  50  km.  The  vollision  partners  of  interest  arc  the  major  constituents  in  the  upper  atmosphere.  The 
recommended  rate  constants  have  been  implemented  in  the  Strategic  High  Altitude  Atmospheric  Radiance 
Code  (SHARC  TO),  which  is  used  to  calculate  background  infrared  radiance. 


20.  DISTRIBUTION/AVAILABILITY  OF  ABSTRACT  21.  ABSTRACT  SECURITY  CLASSIFICATION 

□  unclassified/unumited  0  same  as  rpt  □  dtic  users  Unclassified  _ 


22b.  TELEPHONE  (Include  Area  Code)  22c,  OFFICE  SYMBOL 


22a  NAME  OF  RESPONSIBLE  INDIVIDUAL 


DD  FORM  1473,  84  MAR 


83  APR  edition  may  be  used  until  exhausted.  SECURITY  CLASSIFICATION  OF  THIS  PAGE 

All  other  editions  are  obsolete  UNCLASSIFIED 


AEROSPACE  REPORT  NO. 
TOR-0091  (6078)-2 


RELAXATION  PROCESSES  OF  VIBRATIONALLY  EXCITED  H2O 
IN  THE  MESOSPHERE  AND  THERMOSPHERE 


Prepared 


/  X  •  LcoM___ 


V.  I.  Lang 


Space  and  Environment  Technology  Center 
Technology  Operadons 


Approved 


C.  J.  Rice  O' 

Target  Signatures  and  Backgrounds  Dept. 


/-]  ,  ) 


A.  B.  Christensen,  Principal  Director 
Space  and  Environment  Technology  Center 


’Si- 


1  Aco&aalcu  For 
j  NTT 8  ORi&I 
;  DTiC  TAB 
|  Un/>.'iao"jiiied 

j  Jastiricaciosu 

1 — - ; 

! 

1 


37— _ _ _ 

Distribution/ 

Availability  Codes 
(Avail  and/or 
Dl3t  |  Spsoial 


1 


CONTENTS 


INTRODUCTION . 7 

DATA  REVIEW . 9 

SHARC  RADIANCE  OUTPUT . 15 

SUMMARY . 23 

APPENDIX  A  -  V-R,T  Relaxation  of  H2O  V2  Bend  and  V2  Combination  Levels . 27 

APPENDIX  B  -  V-V  Relaxation  of  H20  by  O2 . 35 

APPENDIX  C  -  V-R,T  Relaxation  of  the  Vj,  V3  Stretching  Reservoir  and 

Combination  Levels . 41 

APPENDIX  D  -  Rate  Constants  for  Spontaneous  Emission  Processes  from 

Vibrationally  Excited  H2O . 47 


FIGURES 


1 .  Vibrational  bands  contributing  to  H2O  infrared  emissions  in  the  upper  atmosphere . 12 

2 .  SH ARC  2.0  limb  calculation  of  H2O  infrared  radiance  using 

old  kinetics  (Ref.  2)  and  revised  constants  (this  review) . 17 

3 .  SHARC  2.0  limb  calculation  of  H2O  infrared  radiance . 18 

4 .  SHARC  2.0  limb  calculation  of  H2O  infrared  radiance . 19 

Al.  V-R,T  relaxation  of  H2O(010) . 28 

B 1 .  V-V  relaxation  of  H20(010)  by  O2 . . . 38 


TABLES 

1.  Rate  Coefficients  for  Relaxation  of  Vibrationally  Excited  H20 . 10 

2 .  Band  Radiance  Comparison . 16 

D 1 .  Einstein  A  Coefficients  for  H2O  Vibrational  Levels . 49 


INTRODUCTION 


This  review  is  the  second  in  a  series  of  reports  intended  to  update  the  state-of- 
know  ledge  regarding  the  temperature  dependence  of  collisional  energy  transfer  processes  by 
potential  radiators  in  the  mesosphere  and  thermosphere.  At  the  altitudes  of  interest  (>  50  km), 
the  radiating  species  are  generally  non-LTE,  that  is,  they  are  not  in  local  thermodynamic  equi¬ 
librium.  For  this  reason,  the  rates  of  spontaneous  emission  processes  that  compete  with  colli 
sional  relaxation  processes  must  be  included  in  any  kinetic  description  of  the  upper  atmo¬ 
sphere.  A  previous  Aerospace  Corporation  report  [TOR-0091  (6078)- 1,  Ref.  lj  described  the 
kinetics  affecting  CO  radiance;  this  report  is  concerned  with  vibrationally  excited  H2O.  A 
discussion  of  the  non-LTE  region  is  contained  in  the  former. 

The  last  major  review  of  atmospheric  kinetics  that  contained  relaxation  rates  for  vibra¬ 
tionally  excited  species  was  that  of  Taylor  in  1974  (Ref.  2).  After  reviewing  the  literature 
through  1991,  we  concluded  that  there  is  still  a  need  for  reliable,  process-specific,  experimental 
rate  data  for  the  collisional  relaxation  of  vibrationally  excited  H2O.  Recent  laser-induced  fluo¬ 
rescence  experiments  have  improved  the  temperature-dependent  data  for  the  self-relaxation  of 
II2O  and  for  the  relaxation  by  0  atoms.  However,  the  temperature  dependences  of  the  relax¬ 
ation  pathways  enabled  by  the  dominant  atmospheric  species,  N2  and  O2,  have  not  been  well 
characterized.  Ambient  temperatures  in  the  mesosphere  and  thermosphere  typieally  range  from 
190  to  1300  K.  We  have  compiled  the  available  data  for  each  relevant  relaxation  partner  in 
order  to  obtain  temperature  functions.  Where  possible,  a  linear  Landau-Teller  form  (Ref.  3) 
was  used.  Where  fewer  data  were  available,  a  simplified  T1/2  temperature  function  was  esti¬ 
mated.  The  results  are  presented  in  Appendices  A  through  C.  A  description  of  the  applicability 
of  the  Landau-Teller  equation  is  given  in  Ref.  1. 

Since  the  rovibrational  line  strengths  for  H2O  have  been  accurately  measured  by  high 
resolution  spectroscopy,  the  vibrational  band  strengths  are  well  established.  The  most  recent 
version  of  the  sum-of- lines  band  strengths,  as  compiled  in  the  HITRAN  1991  listing  (Ref.  4), 
have  been  used  to  calculate  the  spontaneous  emission  rates  from  each  H2O  vibrational  level. 
The  results  are  described  in  Appendix  D. 

The  rate  constants  for  the  energy  transfer  processes  described  here  are  important 
parameters  in  the  calculation  of  radiance  emitted  by  the  upper  atmosphere.  The  values  obtained 
from  this  analysis  have  been  incorporated  into  the  Strategic  High  Altitude  Radiance  Code 
(SIIARCj,  Version  2.0  (Ref.  5).  The  format  for  the  coefficients  of  the  rate  equations  listed  in 
Table  1  has  been  borrowed  from  SHARC.  The  code  was  also  used  to  quantify  the  effect  of  the 
updated  rate  constants  on  the  emissions  from  each  H2O  vibrational  band.  The  results  are 
summarized  in  Table  2  and  illustrated  in  Figs.  2  through  4. 
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DATA  REVIEW 


Some  general  results  and  observations  from  our  review  are  described  in  this  section. 
Specific  references,  data  analyses,  and  rate  constants  are  presented  in  Appendices  A  through  D 
as  well  as  in  Table  1.  A  vibrational  energy  level  diagram  for  H2O  is  given  in  Fig.  1,  depicting 
the  radiative  processes  that  result  in  major  atmospheric  infrared  emissions. 

The  V-R,T  relaxation  of  the  V2  (010)  bending  level  of  H2O  has  been  studied  more  than 
any  of  the  other  collisional  relaxation  processes  relevant  to  atmospheric  H2O.  However,  there 
is  still  an  ambiguity  associated  with  the  effect  of  different  collision  partners  on  the  V2  relaxation 
rate.  The  data  for  the  collision  partners  O2.  N2,  and  Ar  are  combined  in  Fig.  A1  (Appendix  A) 
since  these  are  not  statistically  different  when  the  errors  associated  with  the  measurements  are 
considered.  Previous  versions  of  SHARC  contained  the  same  rate  for  the  relaxation  of  the 
2v2,  3v2,  vi+  v2,  and  vj  +  V3  levels  as  for  the  V2  level.  This  is  not  the  best  approximation  that 
can  be  made  for  all  of  these  higher  energy  levels,  and,  therefore,  we  have  modified  these 
values  as  described  in  Appendix  A. 

The  H2O(010)  level  is  only  39  cm*1  higher  in  energy  than  the  02(1)  level.  This  allows 
for  rapid  V-V  energy  transfer  from  H2O  to  O2  in  a  process  that  is  competitive  with  the  V-R,T 
relaxation  of  1120(010).  The  two  processes  are  difficult  to  separate  experimentally,  but  it  can 
be  accomplished  using  a  carefully  planned  laser-induced  fluorescence  experiment.  However, 
the  only  data  currently  available  for  the  V-V  relaxation  of  the  V2  level  are  acoustic  data,  reported 
with  large  error  bars.  Therefore,  as  shown  in  Fig.  B1  (Appendix  B),  we  have  taken  an 
"average"  of  several  of  these  sets  of  data.  No  temperature  dependence  could  be  established 
over  the  300  600  K  temperature  range  of  the  data.  Our  rate  constant,  based  on  experimental 
data,  is  notably  larger  than  earlier  theoretical  estimates  made  by  Taylor  (Ref.  2). 

Previous  versions  of  the  SHARC  H2O  kinetics  module  already  contained  the  only 
experimentally  measured  room  temperature  rate  constant  for  the  V-R,T  relaxation  of  the  (001) 
and  (100)  reservoir  by  N2  or  O2  (Ref.  6).  Since  there  are  no  newer  data  available  for  other 
temperatures,  we  have  retained  this  constant  for  those  collision  partners.  Current  model 
atmospheres  fe.g.,  1976  Standard,  (Ref.  7)]  generally  indicate  that  the  number  density  of  H2O 
decreases  to  an  insignificant  amount  at  altitudes  where  O(^P)  could  be  a  potential  collision 
partner.  Since  laser  induced  fluorescence  measurements  near  room  temperature  are  available 
for  the  rate  constant  of  V-T  energy  transfer  from  H2CKOOI)  or  (100)  to  0(3P),  these  data  are 
briefly  described  in  Appendix  C  for  use  with  any  non -typical  model  atmospheres  that  may  have 
an  unusually  high  H2O  content  at  high  altitudes. 
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Table  1.  Rate  Coefficients  for  Relaxation  of  Vibrationally  Excited  H2O 

[kj  =  Aj  TP*  exp(>CjT-333  -  E-JT) 
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Figure  1.  Vibrational  bands  contributing  to  H2O  infrared  emissions 
in  the  upper  atmosphere.  Prominent  atmospheric  features  occur 

near  6.4  mm - ,  4.7  |im  — g—  and  2.7  pirn —  —  . 
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Einstein  A  coefficients  for  all  the  relevant  radiative  relaxation  processes  of  vibrationally 
excited  H2O  have  been  calculated  using  the  1991 HITRAN  sum-of  lines  band  strengths,  Sui . 
Changes  to  certain  H2O  rovibrauonal  lines  have  been  made  to  the  HITRAN  listing  since  the 
1986  version  (Ref.  8).  Some  of  the  emission  rates  we  calculated  were  approximately  twice  as 
high  as  those  previously  used  in  SHARC.  Although  the  line  list  for  H2O  has  been  improved 
since  1986,  the  changes  are  not  large  enough  to  cause  such  significant  differences  in  the  band 
strengths.  Rates  calculated  from  the  first  coefficients  of  the  transition  dipole  moments  reported 
in  the  literature  rather  than  from  the  band  strengths  may  have  been  the  source  of  the 
discrepancy.  Dipole  moment  calculations  are  not  straightforward  for  many  H2O  transitions 
since  Fermi  and  Coriolis  interactions  occur  between  closely  spaced  vibrational  levels.  In  most 
cases,  the  full  functional  form  of  the  transition  dipole  necessary  for  calculating  a  band  strength 
is  not  explicitly  given  by  the  authors  who  reported  H2O  line  strengths.  We  have  attempted  to 
standardize  the  emission  rates  using  the  Sui  values  so  that  discrepancies  in  the  application  of  the 
theoretical  parameters  for  the  vibrational  and  rovibrational  interactions  can  be  avoided. 
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SHARC  RADIANCE  OUTPUT 


The  updates  described  above  have  been  incorporated  into  a  new  input  file  for  the  ambi¬ 
ent  H2O  kinetics  module  in  SHARC.  These  rates  are  given  in  Table  1  in  the  format  currently 
used  in  SHARC  (Version  2.0,  Aerospace  Corporation  preliminary  copy).  A  test  calculation 
using  a  SHARC  model  atmosphere  profile  was  performed.  A  subarctic  daytime  model  atmo¬ 
sphere  was  used  since  the  lower  altitude  concentrations  of  H2O  are  larger  than  those  in  the 
Standard  1976  models.  The  results  of  the  test  are  presented  in  Table  2  and  Figs.  2  through  4. 
A  limb  view  at  64  km  tangent  height  was  used  with  an  arbitrary  latitude  and  longitude  of  (0,0) 
chosen  for  illustrative  purposes  only.  Ideally  a  region -specific  equatorial  model  with  large  col¬ 
umn  densities  of  H2O  would  have  been  preferred. 

The  radiances  calculated  by  SHARC  using  the  new  kinetics  scheme  are  generally 
smaller  than  or  equal  to  those  calculated  with  the  old  rates.  The  band  radiances  for  each  H2O 
vibrational  transition  and  the  ratios  (old:new)  are  listed  in  Table  2.  The  decreased  emissions 
are  a  result  of  increasing  the  collisional  relaxation  rates  of  the  upper  vibrational  levels,  or  effec¬ 
tively  quenching  the  excited  molecules  before  they  can  radiate.  This  effect  apparently  domi¬ 
nates  over  the  increases  made  to  some  of  the  spontaneous  emission  rates.  Transitions  originat¬ 
ing  in  the  3v2  level  (centered  at  4667  cm*1, 3072  cm*1,  and  1515  cm*1)  are  a  factor  of  12.6 
times  weaker  after  changing  the  kinetics.  Similarly,  transitions  originating  in  the  2V2  level  are 
5.2  times  weaker.  The  overall  effect  on  the  radiance  in  particular  cnr1  regions  is  illustrated  in 
Figs.  2  through  4.  Our  calculations  indicate  that  the  (020)  -  (000)  band  centered  at  3150  cm*1 
will  be  narrower  due  to  less  intense  radiance  on  the  high  frequency  side  than  previously  pre¬ 
dicted.  This  effect  is  evident  in  Fig.  4.  Such  changes  to  the  more  intense  H2O  infrared  bands 
have  important  consequences  for  atmospheric  sensor  development. 
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Table  2.  Band  Radiance  Comparison 


Transition 

Frequency 

Band  Radiance  (W/sr»cm2)  Radiance  Ratio 

Old  Kinetics  New  Kinetics  01d:New 

(010)-(000) 

1594.750 

0.37277e-6 

0.35466e-6 

1.05 

(020)- (000) 

3151.630 

0.11 184e-8 

0.21412e-9 

5.22 

(100)-(000) 

3657.053 

0.87024e-9 

0.47780e-9 

1.82 

(001)-(000) 

3755.930 

0.35921e-7 

0.36124e-7 

0.994 

(030M000) 

4666.793 

0.1 1070e-l  1 

0.87792e-13 

12.6 

(110)-(000) 

5234.977 

0.63862e-ll 

0.2491  le-11 

2.56 

(Oll)-(OOO) 

5331.269 

0.207 63e-8 

0.12812e-8 

1.62 

(020)- (010) 

1556.880 

0.39202e-7 

0.75069e-8 

5.22 

(100)-(010) 

2062.303 

0.13254e-9 

0.73330e-10 

1.81 

(OOI)-(OIO) 

2161.180 

0.16949e-8 

0.16950e-8 

1.00 

(030)-(010) 

3072.043 

0.17710e-9 

0.14089e-10 

12.6 

(llO)-(OlO) 

3640.227 

0.3901 8e- 10 

0.15291e-10 

2.55 

(Oll)-(OIO) 

3736.519 

0.82490e-8 

0.50967e-8 

1.62 

(030)-(020) 

1515.163 

0.32288e-8 

0.25724e-9 

12.6 

Pure  Rotational  Lines 

0.99151e-10 

0.99151e-10 

1.0 

Total  Overlap  Corrected  Bandpass  Radiance  0.46488e-6 

0.40686e-6 

1.14 

Model  Atmosphere:  DAYSUBAR.DAT  Limb  View,  64  km  Tangent 

Lat.  0°,  Long.  0° 

Solar  Zenith  Angle  0° 
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REVISED  KINETICS 
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Figure  3.  SHAFLC  2l0  limb  calculation  of  H2O  infrared  radiance, 
Main  feature  in  this  region  is  the  (030)  -  (000)  rovibrational  band, 


SUMMARY 


While  there  is  still  room  for  improvement  in  experimental  measurements  of  the 
H2O(010)  vibrational  relaxation  processes,  some  progress  has  been  made  since  1974  when 
several  of  the  rates  reported  by  Taylor  were  based  on  theoretical  inferences. 

Essentially  all  of  the  individual  collisional  relaxation  processes  need  further  experimen¬ 
tal  investigation  as  a  function  of  temperature.  Laser-induced  fluorescence  techniques  would  be 
particularly  useful  in  verifying  previous  acoustic  measurements  as  well  as  accumulating  data  in 
the  temperature  region  below  the  shock-tube  limit  and  above  room  temperature.  The  V-V  and 
V-T  contributions  to  the  overall  relaxation  of  the  V2  bending  level  by  O2  still  need  to  be 
quantified. 

Since  the  vibrational  line  strengths  for  the  H2O  emission  bands  used  in  radiance  codes 
such  as  SHARC  have  been  accurately  measured  by  high-resolution  infrared  spectroscopy,  the 
spontaneous  emission  coefficients  are  well  established. 
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APPENDIX  A 


V-R,T  RELAXATION  OF  H20  v2  BEND  AND  v2  COMBINATION  LEVELS 
a)  Relaxation  of  v2  by  02  N2 

The  two  predominant :  x  ..  the  mesosphere  that  can  act  as  collision  partners  for 
vibrationally  excited  H20  <>  •<„  2  .  ,  The  non-radiative  relaxation  from  the  v2  level  can  be 
expressed  as 

J;  , 

M  +  H20  (010)  +  H20  (000)  M  =  N2, 02  (Al) 

The  H20  kinetics  module  .n  earlier  versions  of  SHARC  contained  the  following  tem¬ 
perature  dependence  for  the  V-T  (or  V-  R,T)  relaxation  of  the  v2  level  of  H20: 

kioId  forward  =  5.37e-10  exp(-69.9i  8/T-333)  200  -  400  K  (A2) 

This  rate  equation  is  from  the  1974  review  article  by  Taylor  (Ref.  Al)  and  has  a  quoted  uncer¬ 
tainty  of  an  order  of  magnitude.  The  original  research  on  which  this  Arrhenius  equation  is 
based  was  that  of  Monk  (Ref.  A2),  who  only  measured  rates  at  295  K,  using  a  sound  absorp¬ 
tion  technique.  It  appears  that  Taylor  derived  an  average  rate  constant  for  the  relaxation  of 
II2O(010)  by  N2  '02  from  two  room  temperature  values  (which  differ  by  a  factor  cf  23)  and 
then  applied  a  theoretical  temperature  dependence.  That  O2  and  N2  are  equally  efficient  part¬ 
ners  is  a  reasonable  assumption  for  the  V-R,T  relaxation  of  H2O(010).  However,  Monk’s  data 
are  too  widespread  to  support  this  assumption  (see  Fig.  Al).  Another  problem  lies  in  the  rela¬ 
tive  magnitudes  of  Monk's  data.  It  is  inherently  difficult  to  separate  the  near-resonant  V-V 
energy  transfer  from  H2O(010)  to  O2(100)  from  the  V-T  relaxation  of  H20(010)  by  02. 
However,  if  the  V-T  contribution  to  the  overall  relaxation  process  is  approximately  the  same  as 
the  V-T  rate  with  N2,  this  fraction  could  be  subtracted  from  the  overall  02  rate.  The  major 
inconsistency  in  Monk's  data  is  the  fact  that  the  V-T  rate  with  02  is  reported  to  be  23  times 
slower  than  that  with  N2 .  We  would  expect  the  former  to  be  larger  by  a  small  amount. 

In  all  of  the  low  temperature  studies  (295  -  600  K)  reported  in  the  literature  (Refs.  A2, 
A3  and  A4),  the  acoustic  techniques  used  to  derive  relaxation  rate  constants  for  Reaction  Al 
depend  on  the  validity  of  approximations  made  for  several  concurrent  energy  transfer  processes 
in  a  mixture  of  gases.  The  experiments  were  generally  designed  to  study  the  dependence  of 
Napier  frequencies  (frequency  of  maximum  sound  absorption  per  meter)  as  a  function  of  O2 
and  H20  concentrations  in  gaseous  mixtures  such  as  air.  While  only  a  few  rates  were  actually 
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C  O 


Figure  Al.  \^-R,T  relaxation  of  H 2O(010).  Data  used  in  present 
:  Kung  and  Center  (Ref.  A6),  N2B,  Ar»;  Bass  et.  al.  (Ref.  A3), 
2  •;  Keeton  and  Bass  (Ref.  A4),  N2  or  Ar  A;  Bass  and  Shields 
(Ref.  A5),02+.  Data  used  bv  Taylor  (Ref.  Al)  and  not  included  in 
present  fit:  Monk  (Ref.  A2),  02  A  ,  N2  □. 


measured,  five  to  seven  reversible  energy  transfer  processes  were  needed  to  stoichiometrically 
balance  the  Napier  problem  and  to  fit  the  experimental  frequencies  to  a  concentration  curve. 
Examination  of  the  procedures  used  by  several  investigators  reveals  that  those  employed  by 
Bass,  Keeton,  and  Williams  (Ref.  A3)  are  more  self-consistent  than  those  of  earlier  investiga¬ 
tors  (Monk,  Ref.  A2).  However,  the  conclusions  drawn  by  Bass  et  al.  are  questionable 
because  they  erroneously  plotted  data  from  other  sources  (e.g.,  Kung  and  Center,  Ref.  A6)  as 
well  as  from  their  own  tables. 

Without  going  into  the  details  of  the  various  approximations  made  by  different  investi¬ 
gators,  we  report  here  the  most  consistent  of  the  results  from  the  acoustic  experiments  (Ref.  A3, 
A4  and  A5)  for  400, 500  and  600  K.  However,  the  individual  data  points  from  these  references 
are  not  necessarily  experimentally  unique.  The  author  Bass  was  involved  in  each  of  these  studies 
and  has  incorporated  the  actual  measurements  into  an  iterative  fitting  procedure  to  compute  the 
temperature  dependence  giver,  in  Ref.  A3.  We  have  combined  the  accumulated  acoustic  data 
with  high  temperature  shock  .uoe  data  for  the  collisional  relaxation  rates  of  H2O  (010).  The  latter 
data  were  taken  from  the,  original  paper  by  Kung  and  Center  (Ref.  A6).  A  statistical  energy 
distribution  using  the  heat  capacity,  Cv,  and  the  vibrational  temperature,  0j,  was  applied  to  the 
experimentally  measured  rates  to  correct  for  bulk  heating  inherent  in  the  shock  tube  method.  The 
authors  estimate  that  their  rate  constants  are  acc  irate  to  within  a  factor  of  2. 

In  summary,  it  is  recommended  that  the  functional  form  of  the  temperature  dependence 
for  the  H2O  (010)  V  T  relaxation  be  updated  to  reflect  the  experimental  shock  tube  data  and  the 
consensus  of  the  acoustic  data.  In  our  analysis,  we  assumed  that  the  V-T  transfer  efficiencies 
of  N2  and  O2  are  approximately  equal.  We  used  the  data  for  the  relaxation  of  H2O  by  Ar  as 
well  as  by  O2  and  N2  since  the  molecular  mass  of  O2  lies  between  that  of  Ar  and  N2.  The 
temperature  dependence  obtained  from  fitting  the  combined  data  to  a  Landau-Teller  function 
(see  Fig.  Al)  is 


kincw  =  1.13  x  10-10  exp(-52.5/T-333)  295  -  2000  K  (A3) 

TheT  333  coefficient  was  obtained  from  the  slope  of  Fig.  Al  while  the  pre-exponential  factor 
is  equal  to  the  exponential  of  the  ordinate  intercept.  At  300  K,  the  rate  constant  from  Taylor’s 
function  is  1.5  x  10" 14  cm3  molecule*1  s*1  compared  with  our  new  value  of  4.4  x  10' 14  cm3 
molecule-1  s  L  However,  this  factor-of-3  difference  decreases  as  the  temperature  increases 
from  room  temperature  up  to  ~1500  K. 

The  rate  for  the  corresponding  reverse  process  (i.e.,  collisional  excitation  of  H2O  in  the 
V2  mode)  includes  the  energy  difference  [E,  =  AE  (cm1  x  1.44  K/cm-1)]  between  ground  state 
H2O  and  H20  (010).  That  is, 
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kincw  reverse  =  1.13  x  10'10  exp(-52.5/T-333  -  2296.4/T) 


(A4) 


where  the  band  V2  band  center  frequency  is  1594.74  cm-1.  Since  the  SHARC  2.0  code 
accesses  the  H2O  States  file  for  the  E(  variable  for  the  V-R,T  relaxation  processes,  these  are  not 
listed  in  Table  1.  The  energies  of  the  vibrational  bands  can  also  be  found  in  the  1991 
HITRAN  (Ref.  A7)  line  list. 

Since  the  low  temperature  rate  constants  for  the  relaxation  of  H2O  by  N2  and  O2  con¬ 
tain  so  much  uncertainty  (~  order  of  magnitude)  it  would  be  worthwhile  to  measure  the  rates 
over  the  180  -1000  K  temperature  range  using  laser-induced  fluorescence  (LIF).  The  bulk 
methods  previously  employed  for  the  study  of  V2  lack  the  capability  of  selective  excitation  of  a 
particular  level.  LIF  has  been  used  successfully  to  study  other  collision  partners  and  other 
vibrationally  excited  levels  of  H2O  (e.g.,  Refs.  A8,  A9). 

b)  V-R,T  Relaxation  of  the  2v2  and  3V2  Levels 

In  previous  versions  of  SHARC,  ki  for  the  V-R,T  relaxation  of  the  H2O  (010)  level  by 
N2  and  O2  was  also  used  for  the  hot  emission  bands,  3V2  -  2V2  and  2V2  -  V2.  Using  the  har¬ 
monic  oscillator  approximation,  the  rates  for  these  V-T  relaxation  mechanisms  of  H2O  would 
be  3  and  2  times  greater,  respectively,  than  that  for  the  \i  -  0  fundamental  process.  The 
revised  rate  expressions  for  collisional  relaxation  with  equal  efficiency  by  M  =  O2  or  N2  are 

M  +  H2O(020)  =  M  +  H2O(010); 
k2ncw  forward  =  2.26  x  10'10  exp(-52.5/T-333)  and 
k2ncw  reverse  =  2.26  x  10'10  exp(-52.5/T-333 . 2241.9/D 

and 

M  +  H20(030)  =  M  +  H20(020); 
k3ncw  forward  =  3-39  x  lO'^  exp(-52.5/T-333)  and 
k3ncw  reverse  =  3.39  x  10**0  exp(-52.5/T-333 . 2181.8/T) 

The  band  center  frequencies  of  1594.74  cm*1,  3151.63  cm-1  and  4666.8  cm'1  for  the  V2, 

2v2,  and  3v2  levels,  respectively,  were  used  to  obtain  the  energy  difference  terms  for  the 
reverse  processes. 

Since  the  rate  coefficients  for  the  overtone  transitions  from  3v2  and  2vi  to  the  ground 
state  are  assumed  to  be  much  smaller  than  those  for  transfer  of  a  single  quantum  of  vibrational 
energy,  they  are  not  included  in  our  overall  kinetics  scheme. 


(A5) 

(A6) 


(A7) 

(A8) 
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Experimental  evidence  for  the  validity  of  using  the  harmonic  oscillator  approximation 
for  the  case  of  H2O  vibrational  relaxation  has  been  found  in  laser- induced  fluorescence  studies 
(Ref.  A8  and  A 10).  This  is  worth  mentioning  here  because  it  might  be  expected  that  strong 
dipole  interactions  or  hydrogen  bonding  by  H2O  would  result  in  non-ideal  gaseous  mixtures. 
The  two  processes  involving  H2O  self-relaxation, 


and 


ka 

H20(020)  +  H20  -  H2O(010)  +  H20 
kb 

H2O(010)  +  H20  -  H20  +  H20, 


(A9) 

(A10) 


were  investigated  experimentally.  It  was  found  that  the  rate  of  the  relaxation  from  2v2  -v2  was 
greater  than  that  from  v2  -  0  and  potentially  twice  as  fast.  The  following  limits  were  assigned, 
2ka>kb>ka  (All) 

The  channel 

H20(020)  +  H20  -  H2O(0 10)  +  H2O(0 1 0)  (A  12) 

was  found  to  be  insignificant  in  comparison.  In  Ref.  A8,  the  relaxation  of  H20  by  rare  gases 
was  not  studied  in  as  much  detail  as  the  self-relaxation  case. 


c)  V-R,T  Relaxation  of  v2  Combination  Levels  (vi+  v2,  v2+  V3) 

For  the  SHARC  reactions  involving  relaxation  of  one  quantum  of  v2  from  a  combina¬ 
tion  level,  the  rate  used  for  relaxation  of  the  v2  level  (kj)  can  be  used.  The  energy  transferred  is 
comparable  in  these  three  cases  (1574  -  1593  cm*1),  and,  therefore,  the  increase  in  rate  as  a 
function  of  temperature  should  also  be  similar. 


M  +  H20(1 10)  -  M  +  H2O(100); 
k5  forward  =  1-13  x  10-10exp(-52.5/T-333) 

(A  13) 

M  +  H2O(100)  -  M  +  H20(110); 

k5  reverse  =  1-13  x  lO1©  exp(-52.5/T-333 . 2272.2/T) 

(A  14) 

M  +  H2O(01 1)  -  M  +  H2O(001); 
k6  forward  =  U3  x  1<H0  exp(-52.5/T-333) 

(A  15) 

M  +  H2O(001)  -  M  +  H2O(01 1); 

k6  reverse  =  U3  x  10-»0  exp(-52.5/T-333 . 2268.5/T) 

(A  16) 
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The  energies  of  the  vi+  V2,  Vj,  V2+V3,  and  V3  band  centers  are  5235.0  cm'1, 3657.08  cm-1, 
5331.279  cm*1,  and  3755.92  cm1,  respectively.  The  band  center  frequencies  were  used  to 
obtain  the  energy  level  discrepancies  for  the  reverse  processes. 
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APPENDIX  B 


V-V  RELAXATION  OF  H20  BY  02 

a)  Relaxation  of  the  H2O(010)  Level 

The  02(1)  level  is  only  38  cm-1  lower  in  energy  than  the  H2O(010)  level.  This  coinci¬ 
dence  facilitates  rapid  resonant  V-V  energy  transfer  from  H2O(010)  to  02.  That  is, 

ki 

H20  (010)  +  O2(0)  =  H20(000)  +  02(1)  (Bl) 

This  process  competes  with  the  relaxation  of  H20  (010)  via  the  V-(R),T  channels  described  in 
Appendix  A.  There  have  been  no  experiments  to  date  designed  to  specifically  separate  the  V-V 
and  V-R,T  contributions  from  the  total  relaxation  rate  of  the  2v2  H20  level  by  02.  This  would 
be  a  difficult  task  even  using  laser-induced  fluorescence  data  since  both  amplitude  subtraction 
as  well  as  decay  time  separation  is  involved.  Alternatively,  if  the  assumption  that  H2O(010) 
relaxation  by  N2  via  a  purely  V-R,T  mechanism  is  approximately  equal  to  the  V-R,T  compo¬ 
nent  of  the  relaxation  by  O2,  a  reasonable  approximation  to  the  latter  rate  could  be  made.  The 
V-R,T  contribution  could  then  be  subtracted  from  the  overall  Cerate,  leaving  the  V-V  rate. 

This  assumption  is  valid  because  there  are  no  resonant  V-V  channels  available  with  N2. 

For  many  kinetic  applications,  the  rate  determining  step  is  the  most  important  process. 
Since  relaxation  via  a  V-V  energy  transfer  is  often  much  faster  than  relaxation  through  a  V-T 
channel,  the  faster  rate  is  usually  neglected.  However,  in  modeling  systems  with  high  vibra¬ 
tional  temperatures,  such  as  H20  in  the  thermosphere,  the  V-V  fundamental  relaxation  rate  of 
H2O(010)  is  used  to  infer  rates  for  four  other  reversible  relaxation  processes  from  higher 
vibrational  levels  of  H20  by  02.  Each  process  involves  the  transfer  of  one  quantum  of  v2 
energy.  Overall,  the  five  processes  involved  could  significantly  affect  the  radiance  from  the 
H20  bands.  For  this  reason,  it  would  be  pref  -able  to  characterize  the  V-V  and  V-T  rate 
constants  explicitly. 

In  Taylor's  1974  review  (Ref.  Bl),  a  temperature-independent  rate  of  1.0  x  10*12  cm3 
molecule'1  s'1  was  assumed  for  the  forward  rate  of  the  reversible  process,  Bl.  This  rate 
constant  was  based  on  a  limited  amount  of  theoretical  information  (Ref.  B2). 

Using  the  Tuesday-Boudart  hypothesis  (Ref.  B3)  for  the  five  reversible  energy  transfer 
processes  possible  in  a  heated  mixture  of  H20  and  02,  Bass  and  Shields  (Ref.  3, 1974),  and 
later  Bass,  Keeton  and  Williams  (Ref.  4, 1976)  found  rate  constants  for  lq  reverse  at  400, 500 
and  600  K.  These  were  not  direct  measurements,  but  rather  the  results  of  an  iterative  proce¬ 
dure  using  experimental  rates  for  three  of  the  five  processes  in  the  mixture.  The  experimental 
values  were  from  acoustic  measurements.  The  data  from  each  of  the  Bass  references  is  plotted 
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in  Fig.  Bl.  The  figure  iso  contains  a  room-temperature  theoretical  prediction  by  Jones  et  al. 
(Ref.  B5).  Since  seve-  of  the  reported  rates  were  for  kj  reverse,  they  were  converted  to  the 
corresponding  forward  rate  constants  using  the  -E/T  term  (where  Ej  =  55.2)  before  being 
plotted. 

There  is  no  temperature  dependence  evident  over  the  very  limited  temperature  range  cov¬ 
ered  by  the  data  in  Fig.  Bl.  The  errors  associated  with  the  data  vary  from  conservative  estimates 
of  50%  to  an  order  of  magnitude;  therefore,  the  data  did  not  warrant  a  full  Landau-Teller  fitting 
procedure.  All  of  the  data  lie  above  or  equal  to  the  temperature-independent  rate  constant  pro¬ 
posed  by  Taylor  (Ref.  Bl)  and  used  in  previous  versions  of  SHARC.  Taylor's  estimated  con¬ 
stant  is  really  a  lower  limit  over  the  ~  200  -  400  K  range.  The  average  of  9  data  points,  including 
one  for  Taylor's  estimate,  yields  a  rate  constant  of  kj  forward  =  2.8  ±  .4  x  10'12  cm3  molecule*1 
s*1.  As  the  temperature  increases,  the  reverse  process,  relaxation  of  C>2(1)  by  H2O,  will  become 
slightly  more  significant.  The  energy  discrepancy  for  the  reverse  process  is  38.3  cm1  so 
kj  reverse  =  2.9  x  10*12  exp  (-  55.8/T).  However,  apart  from  this  shift  in  the  population 
distribution  between  the  H2O(010)  and  02(1)  states,  a  relatively  small  temperature  dr  pendence  is 
expected  for  the  rates  of  such  resonant  energy  transfer  processes. 

b)  V-V  Relaxation  from  Higher  Vibrational  Levels 

In  previous  versions  of  SHARC,  the  same  lower  limit  to  the  rate  of  H2O  (010)  V-V 
relaxation  by  O2  was  used  for  all  the  V-V  relaxation  processes  of  upper  vibrational  levels 
containing  a  V2  mode. 

Applying  the  harmonic  oscillator  approximation  to  the  rate  of  energy  transfer  from  H2O 
to  O2  by  the  V-V  mechanism,  we  can  multiply  the  rate  for  die  H2O  (V2  -  0)  process  by  factors 
of  two  and  three  for  the  (2v2  -  V2)  and  (3V2  -  2V2)  transitions,  respectively.  In  units  of  cm3 
molecules*1  s*1,  the  assigned  rates  are  as  follows: 

H20(020)  +  02(0)  =  H2O(010)  +  02(1); 

k2  forward  =  5-8  X  10*12 

k2  reverse  =  5.8  x  10*12  exp(-1.3AT)  (B2) 

H20(030)  +  02(0)  =  H20(020)  +  02(1); 

k3  forward  =  8-7  x  10*12 

k3  reverse  =  8.7  x  10*12  exp(58.8/T)  (B3) 

Two  other  near-resonant  V-V  energy  transfer  processes  from  higher  combination  levels  of  H2O 
to  the  02(1)  level  are  possible: 


H2O(110)  +  02(0)  =  H20(100)  +  02(1); 

^3  forward  =  2.9  x  10**2 

k3  reverse  =  2.9  x  10'12exp(-31.6/T)  (B4) 

H2O(011)  +  02(0)  =  H20(001)  +  02(1); 

^3  forward  =  2.8  x  10'  ^2 

k3  reverse  =  2.8  x  10'12  exp(-27.8/T)  (B5) 

The  rate  constants  for  Reactions  B4  and  B5  are  approximately  equal  to  ki  since  the  v2 
quantum  number  of  the  upper  H20  levels  is  one,  and  since  only  one  quantum  of  energy  from 
the  V2  mode  is  transferred. 
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Figure  Bl.  V-V  relaxation  of  H20(010)  by  O2.  Data  from  Bass, 
et.  al.  (Ref.  B4),  •  ;  Jones  et.  al.  (Ref.  B5),  A ;  Bass  and  Shields 
(Ref.  B3),  ■  .  Average  from  this  analysis  — .  Constant  from 
Taylor  (Ref.  Bl),  _ 
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APPENDIX  C 


V-R,T  RELAXATION  OF  THE  vi,  v3  STRETCHING 
RESERVOIR  AND  COMBINATION  LEVELS 

a)  H20(100),(001)  Relaxation  by  O2,  N2 

The  energy  difference  between  the  band  centers  of  the  Vi  and  V3  stretching  levels  of 
H2O  is  99  cm1.  Intramolecular  energy  transfer  between  these  levels  is  facilitated  by  collision 
with  a  third  body, 

kl 

M+  H20(001)  =  H2O(100)  +  M  +  99  cm'1  (Cl) 

Single  collisions  with  a  third  body  are  probably  enough  to  cause  this  intramolecular  energy 
transfer  in  H20  to  occur  at  room  temperature.  The  rapid  equilibration  between  the  (001)  and 
(100)  levels  has  been  observed  by  Finzi  et  al.  (Ref.  Cl).  We  estimate  that  the  rate  constant  kj 
is  gas  kinetic  at  room  temperature.  Based  on  the  collision  cross-sections  of  H20  and  N2,  the 
temperature  dependence  is  estimated  to  be 


k!  =  1.2  x  10’11  T1/2  (C2) 

The  bidirectional  interconversion,  reaction  Cl,  is  rapid  compared  to  that  of  energy  transfer  pro¬ 
cesses  to  lower  vibrational  levels.  Furthermore,  V-T  energy  transfer  from  either  the  Vi  or  V3 
level  to  the  2v2  level  (505  cm'1  or  604  cm'1  lower  in  energy,  respectively)  occurs  at 
essentially  the  same  rate. 

One  rate  constant  is  assigned  to  both 

k3 

H2O(100)  +  M  =  M  +  H20(020)  (C3) 

U 

and  H2O(001)  +  M  =  M  +  H20(020)  (C4) 

Finzi  et  al.  (Ref.  Bl)  measured  laser-induced  fluorescence  from  the  (vi,  v3)  reservoir  to  obtain 
a  room  temperature  relaxation  rate  constant  to  the  H20(020)  level.  The  experimentally  mea¬ 
sured  quantity  is  actually  a  weighted  average  of  k3  and  Lf.  Their  values  for  k3  forward  were  4.6 
d*  1.2  x  lO’1^  cm^  molecule’1  s’1  for  M  =  N2  and  3.3  ±  1.2  x  lO'1^  cm^  molecule'1  s’1  for  M 
=  C>2-  Efficiency  coefficients  of  1.0  and  0.7  are  used  to  differentiate  between  the  N2  and  02 
rate  constants  in  Table  1  and  in  SHARC.  Unfortunately  no  information  regarding  the  tempera- 
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ture  dependence  of  these  rate  constants  is  yet  available.  Kung  and  Center  (Ref.  C2),  who 
investigated  other  HiO  vibrational  levels  in  their  shock  tube  experiments,  mention  an  overall 
rate  for  the  relaxation  at  the  stretch  reservoir,  but  they  included  processes  whereby  2V2  and  V2 
were  relaxed  to  the  ground  state.  They  made  the  assumption  that  the  latter  step  was  rate 
determining;  this  was  later  disproven  in  the  work  of  Finzi  et  al.  Earlier  rate  constants  sug¬ 
gested  by  Taylor  (Ref.  C3)  or  Sarjeant  (Ref.  C4)  for  the  (vi,  V3)  reservoir  relaxation  were  not 
based  on  experimental  evidence  specific  to  either  Reaction  C2  or  C3. 

The  Finzi  et  al.  room  temperature  rate  constants  of  4.6  x  10'13  and  3.3  x  10"13  cm3 
molecule  1  s1  for  N2  and  O2.  respectively,  were  already  incorporated  into  the  kinetics  scheme 
of  earlier  versions  of  SHARC,  and,  therefore,  we  have  not  recommended  any  changes. 

b)  Relaxation  of  H2O(110)  and  H2CKOII)  Combination  Levels  by  O2  and  N2 
The  following  V-R,T  energy  transfer  processes  from  combination  levels  are  included  in 
the  H2O  kinetics  scheme  used  in  SHARC: 

k5 


M+  H20(1 10)  =  M  +  H20(030) 

(C5) 

and 

k6 

M+  H2O(01 1)  =  M  +  H20(030) 

(C6) 

with  M=  O2  or  N2.  The  energy  difference  between  the  V]+  V2  level  and  the  V2+  V3  level  is 
only  96  cm1.  Therefore,  rapid  energy  exchange  will  occur  between  the  two  levels,  and  they 
can  be  treated  as  a  reservoir.  Furthermore,  processes  4  and  5  are  analogous  to  processes  2  and 
3  described  in  Part  (a)  of  this  appendix,  that  is,  k2  =  k3  =  ks  =  k^.  A  common  room  tempera¬ 
ture  rate  constant  is  recommended.  Again,  the  rate  constants  for  N2  and  O2  are  differentiated 
by  the  efficiency  factors  of  1.0  and  0.7. 

c)  V-T  Relaxation  of  the  vi  and  V2  Levels  by  Atomic  Oxygen 

The  rate  constant  for  processes  C3  and  C4  with  0(3P)  as  the  collision  partner  was  mea¬ 
sured  by  Zittel  et  al.  (Ref.  C5)  at  315  K  and  520  K  using  laser-induced  fluorescence.  The 
experimental  values  were  3(±  2)  x  10'12  cm3  molecule'1  s'1  and  5(±  3)  x  10'12 
cm3  molecule  1  s  *,  respectively.  There  is  insufficient  data  to  establish  a  tempuature  trend  so 
an  average  rate  constant  of  4  x  10'12  cm3  molecule1  s'1  is  recommended.  In  Table  1,  the  rate 
constant  is  reported  as  a  collision  efficiency  parameter  to  the  rate  constant  for  Reaction  Cl. 

The  1976  Standard  model  atmosphere  (Ref.  C6)  indicates  that  at  altitudes  where  0(3P) 
becomes  a  significant  atmospheric  constituent,  the  number  density  of  H2O  has  decreased  to  the 
extent  that  collisional  relaxation  of  vibrationally  excited  H2O  by  atomic  oxygen  is  not 
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competitive  with  similar  mechanisms  involving  N2  or  O2.  This  conclusion  is  based  on  combin¬ 
ing  typical  number  densities  for  H2O  ana  O(^P)  with  the  average  rate  constant  for  processes  C3 
and  C4. 
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APPENDIX  D 


RATE  CONSTANTS  FOR  SPONTANEOUS  EMISSION  PROCESSES 
FROM  VIBRATIONALLY  EXCITED  H20 

The  1991 HITRAN  sum-of-lines  band  strengths  (Ref.  Dl)  given  in  Column  3  of  Table 
D1  were  used  to  calculate  the  Einstein  A  coefficients  for  spontaneous  emission  (in  units  of  s'1). 
The  formula  used  for  the  calculation  is  that  used  in  the  SHARC  2.0  manual  (Ref.  D2)  and  in 
Ref.  D3.  That  is, 

a  _  8tov2Sui  exp[Et/205.727] 
gu  l-exp[-v/205.727] 

where  v  is  the  transition  frequency  in  cm1  and  Ej  is  the  lower  state  energy  in  cm1.  These  and 
the  other  variables  are  described  in  references  D2  and  D3  as  well  as  in  the  first  report  of  this 
series  (TR-0091(6078)-1).  The  units  for  Siu  are  cm'Vmolecule/cm2  or  cm/molecule. 

A  comparison  of  our  calculated  A  values  and  those  used  in  earlier  versions  of  SHARC 
indicated  several  discrepancies.  Several  of  these  differences  were  <  10%.  However,  the 
vi  -  0,  3v2  -  0,  2V2  -  V2,  vi  -  V2,  and  V3  -  v2  Einstein  A  values  were  in  much  worse  agree¬ 
ment.  It  appears  that  some  of  the  spontaneous  emission  rate  constants  used  in  SHARC  could 
have  been  derived  frcm  ir.e  dipole  moments  reported  in  the  o.iginal  references.  The  Einstein  A 
coefficient  calculated  .rorr.  the  basic  formula  containing  the  transition  moment  Rnm  and  the  fre¬ 
quency  of  the  transition,  v,  is  (Ref.  D4) 

A=  — 64zi£Yl_jRnm|2  (D2) 

(47t£0)  3 he3 

The  latter  calculation  is  complicated  when  there  are  Termi  interacdons  or  Coriolis  effects 
occurring.  This  is  the  case  for  the  two  sets  of  triad  levels  in  H20.  The  interacting  levels  are  vi, 
V3,  and  2v2  in  one  set;  and  3V2,  V1+V2,  and  V2+  V3  for  the  second  set.  Therefore,  it  is  neces¬ 
sary  to  use  a  transformed  dipole  moment  operator  that  includes  a  rotational  operator,  and  a  set 
of  vibration -rotation  wavefunctions  for  all  three  levels  of  the  triad.  The  approximation  that 
IRnm|2  =  |p|2  results  in  the  A  coefficients  in  Column  6  of  Table  Dl.  We  did  not  pursue  the 
detailed  calculations  involved  for  some  of  the  other  vibrational  transitions  via  the  dipole 
moment  method. 
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Since  Rothman's  band  strengths  contain  the  best  literature  values  to  date  for  the 
rovibrational  line  strengths  and,  thus,  the  sum-of-line  strengths,  our  recommendation  is  to 
standardize  the  Einstein  A  coefficients  in  the  code  by  calculating  them  from  Eqn.  Dl. 
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Table  Dl.  Einstein  A  Coefficients  for  H2O  Vibrational  Levels 
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Superscripts  |e.g.,  (DB)J  refer  to  References  Dl  through  Dll. 
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